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Abstract 

We present a BFKL phenomenology of structure functions in the framework of the gener- 
ahzed BFKL equation for the dipole cross section proposed by us recently. We show that 
the HERA measurements of excitation of charm at <10 GeV^ allow determination of 
the intercept of the BFKL pomeron. A comparison with the experimental data suggests 
a large contribution from the nonperturbative cross section for colour dipoles of large size r. 
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The ^'-matrix of diffractive scattering is exactly diagonalized in terms of the dipole cross 
section [1-3] which satisfies the generahzed BFKL (Bahtskii-Fadin-Kuraev-Lipatov) equation 
derived in [4-7]. In [6,7] we have determined the intercept of the pomeron Aip and found the 
pomeron dipole cross section in the realistic model with the running QCD coupling as{r) and 
the finite correlation radius Rc for perturbative gluons. With the fic = ^/ Rc = 0.75GeV and 
the infrared freezing of the strong coupling as{r > Rf) — a^J^^ — 0.8 we found Ajp = 0.4 
[6]. 

In this paper we discuss the coresponding BFKL phenomenology of deep inelastic scat- 
tering at large ^, where x is the Bjorken variable. The starting point of our analysis is the 
dipole-cross section representation for the virtual photoabsorption cross section [1] 

c7t,l{^*N,^,Q') = £dz J d'f\^TA^,r)\M^,r), (1) 

where r) is the dipole cross section for the colour dipole of size r, is the virtuality 
of the photon, ^ = log(^) and xq ~0. 1-0.01 corresponds to the onset of the leading-log(^) 
approximation. In the perturbative part of the dipole cross section a^^\^,r), which is a 
solution of our generalized BFKL equation, the gluon correlation radius Rc emerges as the 
principal parameter which controls its r dependence and the absolute normalization. At 
r > Rc the dipole cross section contains also the nonperturbative component, which is of 
major concern from the point of view of the phenomenology to be presented in this paper. 
The wave functions of the (T) transverse and (L) longitudinal virtual photon were derived 
in [1] and read 

|*T(.,r)p = pj|E4{[-' + (l-^)Vi^iM^ + mJXo(£rn , (2) 

\^L{z,r)\^ = g|E4e? Q'Al - zfK,{erf , (3) 

where Ki{x) is the modified Bessel function, e/ is the quark charge in units of the electron 
charge, = z{l — z)Q'^ + rrij, rrij is the quark mass and z is the fraction of photon's 
(lightcone) momentum carried by the quark of the qq pair (0 < z < 1). Then, the sea 
structure function is calculated as F2{x,Q'^) — Q'^ [ctt + ctl] /{An'^aem}- The dipole cross 
section a{^, r) is universal, the and the flavour dependence of structure functions only 
comes from I P . 



The derivation [1] of (|1|) uses the observation that at ^ ^ 1, the virtual photoabsorption 
can be viewed as interaction with the target proton of the multipartonic hghtcone Fock 
states (gg, qqg---) of the photon, which are formed at large distance Az ~ upstream 
the target. In [4] we gave the regular procedure for construction of the corresponding 
multiparton lightcone wave function and of the multiparton cross section. The Fock states 
with n soft gluons give the oc contribution to the total photoabsorption cross section, 
which can be reabsorbed into the energy dependent perturbative dipole cross section 

^^^*He,r) = E3^nWr, (4) 

n=0 

which satisfies the generalized BFKL equation [2-5] 

?^^^K..a'"\tr), (5) 
where in terms of the expansion (H) the kernel JC is defined by 



(^n+i{r) = /C (g)a„(r) 



.2 

gs{Ri)Ki{nGpi)— - gs{R2)Ki{fiGp2) — 

Pi P2 



.OTniPl) + <yn{p2) - O-n(r)] . (6) 



Here Rc = ^/pg is the correlation radius for perturbative gluons, P2 = Pi — r, Ri = 
min{r. Pi}, gs{r) is the effective colour charge, 

where (3q = 11 — ^Nf = 9 for Nf = 3 active flavours, C ~ 1.5 [1] and we impose the infrared 
freezing as{r > Rj) = al/'^'* = 0.8 [6]. The correspondance to the original BFKL equation 
[8] was discussed in detail in [4-7]. 

Once the dipole cross section @) is known, Eq. (|l]) enables the parameter-free calculation 
of the proton structure function at small x, which is the subject of the present communi- 
cation. Below we discuss the salient features of the emerging BFKL phenomenology of 
deep inelastic scattering at small x. We start with the discussion of the nonperturbative 
component of the dipole cross section, which dominates at large r and at moderately large 
-. Then, we present calculations of the proton structure function in the kinematical range 



of HERA and find good agreement with the data from ZEUS and HI experiments [9,10]. 
We comment on determination of the pomeron intercept from the charm structure function. 
We also give an estimate of the gluon distributions which correspond to these structure 
functions and give predictions for R = ai/cTT- 

We start presentation of our resuhs with the description of the colour dipole cross secton. 
In the perturbative QCD, the growth of the dipole cross section comes from the logarithmic 
growth of the multiplicity ng of perturbative guons in the photon generated from the parent 
qq Fock state Ug oc J^°dzg/zg = Equations (||,|) describe the perturbative part of the 
dipole cross section: the kernel /C is directly proportional to the probability of radiation 
of the perturbative gluons [4-7], and Aag = |[cr(pi) + a{p2) — cr(r)] shows how much the 
interaction cross section is changed for the presense of the perturbative gluon (here pi^2 is 
the separation of the gluon from the (anti) quark, p2 = pi — r and a{r) is the dipole cross 
section for the parent gg-dipole of size r). To higher orders, the gluons themselves become 
the source of the subsequent generation of gluons, so that Fock states with n soft gluons give 
oc contribution to the dipole cross section. As a boundary conditon at a;o = 3 ■ 10~^, one 
can start with the two-gluon exchange dipole cross section for interaction with the nucleon 
target [1] 

Mr) = ^ / -^.-^^as{e)as{^')[l - 0,(3^)] 1 - exp(.A:f) , (8) 

where Gp{q^) is the charge form factor of the proton, = max{A;^, ^}. 

The rightmost singularity in the complex j-plane corresponds to the pomeron cross 
section crip(,^, r) = ajp{r) exp(A]p^) . The crucial point of the BFKL phenomenology is that 
the pomeron cross section o"ip(r) is not, and can not, be confined to the small-r region, which 
is due to the so-called diffusion property of the Green's function of the BFKL equation [8]. 
The underlying physics is very simple: even when one starts with the small-size beam and 
target colour dipoles, radiated gluons stick out of the parent dipoles the transverse distance 
~ Rc apart, and it is interaction of these gluons which controlls the total cross section. For 
instance, the pomeron intercept Ajp is predominatly controlled by the behavior of (Tip(r) in 
the semi-perturbative region of r ~ i?c [6]. Consequently, the contribution from the infrared 
region of large r, and the infrared regularization of the perturbative QCD, are at the heart 



of the BFKL phenomenology. 

Motivated by the lattice studies of the correlation function of perturbative gluons [11], we 
take hg = 0.75GeV. With such a large value of fie, the driving term (|]) of the perturbative 
dipole cross section cr''^^\^ = 0, r) = cro(r) flattens at large r at ~ 8mb (Fig. 1). One 
can easily anticipate the nonperturbative mechanism of interaction of large, r > Rc, colour 
dipoles. Here we make an assumption that the perturbative and nonperturbative cross 
sections are additive [12] 

a(e,r) = (T(P*)(e,r) + (T("P*)(r). (9) 

Because the rise 

ofthea(P*)(^, 

r) is due to the rising multiplicity of the perturbative gluons in 
the photon wave function, we make a plausible assumption that the nonperturbative dipole 
cross section does not rise with energy. We also assume the dominance of the perturbative 
cross section, a^"-^^y a^P^\r) ^ 1 at r <^ Rc- In Fig. 1 we show our choice of a^^^*\r). 

Different processes probe cr(^,r) at different r, which allows to disentangle the pertur- 
bative and nonperturbative dipole cross sections experimentally. For instance, we shall see 
that the charm production predominantly probes a^^*\$,,r). The above choice of a("P*)(r) is 
already strongly constrained by various experimental data. As a matter of fact, we require 
that the sum of the perturbative and nonperturbative cross sections (t{^ = 0, r) roughly 
reproduces the colour dipole cross section of Ref. 1, which was constrained to reproduce 
atot{TTN) ^ 25mb. The cross section of Ref. 1 has already been succesfuUy used in the 
calculation of the structure functions at moderate (see also below) and of the nuclear 
shadowing [1,3], of the total photoproduction cross section and of the rate of the diffraction 
dissociation of real and virtual photons [2]. The E665 data on the exclusive leptoproduction 
of the p°-mesons on nuclei [13] and the NMC data on the J/\E' photoproduction on nuclei 
[14] imply that a{C, = 0, r) decreases from ~ 30mb at r ~(l-1.5)f to ~ 6mb at r ~ 0.4f (for 
the detailed discussion see [15,16]). 

Because the ratio a{^, r)/r'^ only slowly increases towards small r, it is useful to rewrite 
Eq. (|l]) in the form 

F2{x,Q' = - —■ r') . 10 

The divided by contributions of light and charmed quarks to the kernel $r(Q^,r^) = 
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(n'^/Aaem) Jo dz Q'^r'^\'^T,Liz,r)\'^ are shown in Fig. 2. They develope a plateau at large 
Q^, the emergence of which signals the onset of the GLDAP (Gribov-Lipatov-Dokshitzer- 
Altarelli-Parisi) [17] evolution of parton densities. The Fig. 2 in conjunction with Fig. 1 
shows which region of r is probed in deep inelastic scattering at the virtuality Q^. Following 
[1-3], we take the quark masses rriu^d = ISOMeV, m,, = 300MeV, = 1.5GeV, = 
4.5GeV. The structure functions are insensitive to cr(^, r) at r > 2f. Notice, that the 
contribution from r ~ If is the perfectly scaling part of F2{x,Q'^) which persists at all the 
values of [1]. Numerically, we find that the contribution from 

changes form Ft^'\x,Q^ = O.TSGeV^) ^ 0.1 to Fi"^*^(x,g2 = iSGeV^) ^ 0.13. In Fig. 3a 
we show the total sea structure function, in Fig. 3b we show its perturbative part. Because 
of expansion of the plateau in Fig. 2 with rising Q^, the contribution to F2{x,Q'^) from 
^{pt)(^^^^^ takes over at large Q^. The values of structure function already at ^(2- 
4) GeV^ are insensitive to masses of light quarks. Contribution from cr^^*)(^,r) also takes 
over at large ^, but at a; ~ 10~^ the effect of Q^) is quite substantial. 

In Fig. 4 we compare our predictions with the NMC [18], ZEUS [10] and HI [11] structure 
functions. The overall agreement is good; the slight departure from the NMC data points 
at large x is due to the the valence structure function not considered here. Besides, at 
large x the nonperturbative cross section can contain the Regge- decreasing component oc 
exp(— 1^) = \Jx/xq, for the magnitude of which we do not have any educated guess [19]. 
It is somewhat more accurate to use in Eq. + log2;,r), but for the purposes of the 

qualitative phenomenology of the present paper, this rescaling of the energy variable can be 
reabsorbed into renormalization of the boundary condition. 

In [7] we have shown that the energy dependence of a^^^\^,r) at r ~ 0.5Rc follows the 
asymptotic law oc exp(A]p^) starting already at moderate energy. At r <^ the effective 
intercept is substantially larger than Aip beacuse of the double-leading-logarithm effects 
[5,7]. Eq. (15) and Fig. 2 show that the total structure function F2{x,Q'^) receives con- 
tributions from a broad range of r. Besides, Figs. 3a and 3b show that the x-dependence 
of F2{x,Q'^) is strongly affected by the contribution from the nonperturbative cross sec- 
tion, which does not allow the model-independent determination of the pomeron intercept 
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Aip form the total structure fucntion measurements in the kinematical range of HERA. 
However, the same Fig. 2 shows that one can zoom at the perturbative cross section at 
r ~0.15f measuring excitation of charm at ^ 10^, so that in this case the effective in- 
tercept Aeff{x,Q^) = -d\ogF2{x,Q'^)/d\ogx Ajp = 0.4 . At larger Q^, Aeff{x,Q^) 
overestimates Ajp because of the double-leading-logarithm efects [7]. Our predictions for 
the charm structure function are shown in Fig. 3c (The threshold effects in excitation of 
heavy flavours at small [20] can approximately be taken into account by the rescaling 

2 

$, =^ { -I- log(l -I- ^) ) . Once the pomeron intercept and the perturbative dipole cross section 
are determined, one can separate the nonperturbative contribution from the high-accuracy 
data on F2{x, Q"^) at relatively low Q'^ ^OGeV^. To this end, the important point is that 
at ^ (10-20) GeV^ the BFKL evolution and the conventional GLDAP (Gribov-Lipatov- 
Dokshitzer-Altarelli-Parisi) [17] evolution complemented by the proper boundary condition, 
yield essentially identical predictions for the F2^*'' {x, Q"^) in the HERA region (the deep 
reasons for this similarity are discussed in detail in [7]). 

Above we only presented the results for the bare pomeron. At large ^ the unitarization 
effects become important. They were discussed to great detail in [4,21]. The unitarity 
bound can best be formulated in terms of the profile function (the partial wave amplitude) 
r(^, r, h) for the dipole cross section 

52 1 

25^J ' ^^^^ 

where h is the impact parameter. The diffraction slope So(^, r) is a slow function of ^ and 
r compared to the rapidly changing (T(^,r). Evidently, the s-channel unitarity constraint 
^{'ii^i^) < 1 will first be violated at large r. To a good approximation, the unitarization 
effects reduce to a (very weakly Q'^ dependent) renormalization factor (1 — 5^) ~ exp(— 5^) 
with 5 ~0. 03-0. 05 [4,21], which is neghgibly small compared to Ajp = 0.4. Furthermore, 
in [4] we have shown that the unitarized structure function still satisfies the hnear GLDAP 
evolution. Consequently, for the purposes of the present crude phenomenology, the unita- 
rization effects can be reabsorbed into the (as yet unknown) intercept of the pomeron. 

The ratio R — (7l/(7t was first evaluated in the framework of the dipole cross section 
technique in [1]. In [21] we have shown that the unitarization effects have little impact on 



<yL/<yT- Our predictions for R for the sea component of the structure function are shown in 
Fig. 6. We find a sort of crossover at a; ~ (1 — 2) ■ 10"'^. At small x, R rises a little with 
Q^, and is essentially flat at > 20GeV^. The knoweledge of R at small is important 
for certain applications like an accurate estimate of the radiative corrections. In Fig. 6b we 
separately show our predictions for R ai < 0.75 GeV^. Notice that \^ oc Q^, and at 
very small we have R oc just by virtue of gauge invariance. 

The relationship between our dipole-cross section representation for deep inelastic scat- 
tering and the more familar parton model (recall that the latter comes along with the 
Weizsacker- Williams reinterpretation of the photoabsorption cross section [17]) is estab- 
lished by [3,4,21,22] 

a(i,r) « -a.(r) / ^ ^ • « ^rV.,,(DG,,(^. ~ -j) . (12) 

where Gi^^.Q"^) = xgi^XjQ"^) is the gluon structure function at a; = xoexp(— ^). Then, 
assuming Nf = 4 active flavours, the differentiation of Eq. (p!OD gives (see also [23]) 

where Cj is the quark charge in units of the electron charge. Strictly speaking, Eq. ([T2|) was 
derived for the perturbative dipole cross section, but one may use Eq. (0) in conjunction 
with Eq. (0) as a definition of the small-x gluon stucture function beyond the perturba- 
tion theory. If we apply the procedure ([T3|) to our BFKL structure functions, then we find 
G{x, Q^) shown in Fig. 7 in comparison with the recent HI [24] and NMC [18] determina- 
tions of G{x, Q^). We have a good overall agreement with the experiment. 

Discussion of the results and conclusions: 

We presented the BFKL phenomenology of deep inelastic scattering at large ^ in the dipole- 
cross section representation. The nover features of our approach, which make it different 
from the related work by other authors [25], are as follows: Firstly, its irrefutable advantage 
is an exact diagonalization of the S-matrix in the dipole-cross section representation. Sec- 
ondly, the general idea, formulated and exploited in our early works [1-3,20-22], was to use 
the universality property of the dipole cross section to fix the normalization and the shape 



8 



of (t({, r) at moderately large energies. Thirdly, the the nonperturbative dipole cross section 
which emerges as an integral part of the formalism, is well constrained by the hadronic cross 
sections [1], the small-Q^ structure functions [1,3], the diffraction dissociation of photons 
[2], the nuclear shadowing [3] and the leptoproduction of vector mesons [15,16]. Fourthly, 
the discussion in [25] centered on matching the GLDAP solutions with the fixed-a^ solu- 
tions of the scaling BFKL equation [8] . Solutions of our generalized BFKL equation, which 
consistently treats the running coupling CKs{r), are dramatically different from the fixed-ag 
solutions, (for the more detailed discussion see [7]). Finally, our infrared regularization is 
consistent with gauge invariance [4-6] 

In our earlier work [1-3,20-22] we tried to describe the whole dipole cross section a{^, r) 
purely perturbatively. As far as the absolute value of cr(^ = 0, r) is concerned, that can be 
achieved at the expense of small fic- Such a perturbative glue always starts with g{x, Q^) oc 
I and always generates the sea q{x,Q^) oc ^, providing a nice continuity from the real 
photoproduction to deep inelastic scattering. The trouble with small /iq is a too rapid a 
buildup of the perturbative glue, which can be tamed, and a good quantitative description 
of the ZEUS and HI structure functions can be obtained, at the expense of starting the 
GLDAP evolution of glue at Ql IGeV^ 4//^ [21] (the unitraization effects were also 
found important to bring the predictons [21] closer to the ZEUS and HI data). There 
is much correspondance between this large value of Qq of Ref. 21 and the small gluon 
correlation radius Rc used in the present paper. The present analysis shows that in order to 
understand the HERA data, one really needs a large nonperturbative dipole cross besides 
the perturbative BFKL cross section. 

We conclude listing the uncertanties and challenges: The pomeron intercept and the 
pomeron cross section are evidently sensitive to the specific choice of the gluon correlation 
radius Rc and of the freezing coupling, and also can change if more sophisticated infrared 
regularization is introduced. Another uncertainty is the choice of xq. Hovever, none of the 
above will drastically modify the salient features of the above presented dipole-cross section 
phenomenology . The further progress in the theory requires the experimental determination 
of the pomeron intercept Ajp, which the HERA experiments on excitation of charm can do. 



9 



and better understanding of the nonperturbative cross section, which also can be probed in 
the HERA experiments. From the theoretical side, more work on the unitarization effects 
is in order. 

Acknowledgements: B.G.Z. is grateful to J.Speth for the hospitahty at IKP, KFA Jiilich, 
where this work was initiated. 
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Figure captions: 

Fig.l - The nonperturbative and perturbative components of the dipole cross section. 

Fig.2 - The kernel ^T{Q'^,r'^) for (the top box) the hght {u,d) and (the bottom box) the 
charmed quarks. 

Fig. 3 - The predicted BFKL structure functions at small x: (a) the total structure function, 
(b) the perturbative part of the structure function, (c) the charm structure function. 
The curves correspond to (from bottom to top) Q^=0.75, 1.25, 1.75, 2.5,3.5,4.5,8.5,- 
15, 30,60,120,240,480 GeV^, respectively. 

Fig. 4 - Comparison of the predicted BFKL structure functions with the data from the ZEUS 
[10], HI [11] and NMC [18] experiments. 

Fig. 5 - Predictions for R = (Xl/ctt- The lower box shows the low-Q^ behavior of R. 

Fig. 6 - Comparison of the gluon structure function evaluated from the predicted BFKL struc- 
ture functions using Eq. (|13]) with the HI [24] and NMC [18] determinations. 
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